Comment on the paper "Bound States in the One-dimensional Hubbard Model" 

F. H. L. Efiler 

Department of Physics, Theoretical Physics, Oxford University, 1 Kehle Road, Oxford 0X1 NP, UK 

F. Gohmann, V. E. Korepin 
Institute for Theoretical Physics, State University of New York at Stony Brook, Stony Brook, NY 11794-3840, USA 



We comment on the preprint [cond-mat/OSOSlOS by D. Braak and N. Andrei We point 



out that the "new" Bethe Ansatz equations presented in are identical to the Bethe equations 
for strings introduced by M. Takahashi for the description of thermodynamics in 1972 Some 
physics suggested in [|lj is incorrect. In particular, all former conclusions made on the basis of the 
string Bethe equations remain valid. 

(I) First let us correct a typo in . The bare S'-matrix for the scattering of an unbound particle with an m-complex 
(formula (17) in should be 

^«(m) _ sin k - (/)[)'"^ - miu/A 
'='^0"' sinfc-0[,"^ +miu/4 

(see eq. (6.59) of |^). Let us also mention that by convention the S'-matrix in |]l| gets inverted upon horizontal 
transposition of the indices. A straightforward calculation shows that the "new" Bethe equations of |jl|| ((BAG), 
eqs. (20)-(22) in 0]) coincide with the standard string Bethe equations [§. This calculation involves adjustment of 
notations and making explicit the so-called A-strings, which are hidden in the notation of [Q (see appendix A). [The 
A strings are configurations of the spin-rapidities involving complex . The A-strings have to be distinguished 
from the so-called fc-A-strings, which are the subject of criticism in [|] (fc-A-strings involve spin-rapidities A-,, as well 
as charge rapidities kj).] Since the "new" Bethe equations in ||l[ coincide with the string Bethe equations, they can 
not lead to new physics. It is no surprise, in particular, that the counting of states mentioned in [Q gives the correct 
number 4^, since this counting is actually identical to the calculation of The Algebraic Bethe Ansatz expression 
at the start of page 2 of Q is incorrect: the Hamiltonian does not act in this way but needs to be replaced by the 
inhomogeneous transfer matrix of the spin problem. 

(II) Let us now demonstrate that the "new physics" discussed at the end of [1] is partially incorrect and partially 
known. Let us start with an incorrect statement of |]l|: in the last part the holon-antiholon excitation for the half- filled 
repulsive Hubbard model is considered. This excitation is known in the literature as "charge singlet" (we will follow 
the discussion given in |^,^). It is described by two holes in the fc-Fermi sea {k^ and k!^) and one k-A string. The 
spin parameter of the fc-A string is denoted by (j){q) in Q and by A' in ||,^. The authors of |]l| state that the relation 

2A' = sin k'^ + sin fc^ (2) 

is not valid (remember that A' = (j^il))- This statement is incorrect: It was shown in that (^) follows from the 
string Bethe equations (which are identical to (20)-(22) alias (BAG) in |^). Independent confirmation can be found 
in 0. We repeat the derivation of (§) in detail in appendix A below: Eq. (||) holds at half-filling! This confirms the 
expression for the S'-matrix obtained in H,^. 

Let us comment on the "new" excitation below half-filling, which is discussed in the last part of . Our comments 
are: (i) the existence of this excitation was well known ||^ and (ii) is of no particular physical significance, since it (a) 
has a gap and (b) is only one of an infinite number of gapful excitations below half-filling. It is correct that eq. (g) 
is not satisfied below half-filling and that as a result the fc-A string is an independent excitation. However the same 
holds for the infinite number of longer fc-A strings and the one-particle excitations! These excitations also exist at 
zero density (zero filling) |^ . It is the generic situation for Bethe Ansatz solvable models that the structure of the 
excitations persists, when going from a trivial to a non-trivial ground state by changing the density (more generally, 
an order parameter). This means in case of the Hubbard model, that all strings as well as the one-particle excitations 
are present at finite density below half-filling. Only their dispersion curves get modified. Similarly, for the XXX spin 
chain the string excitations over the ferromagnetic ground state survive, when the magnetization is diminished. Only 
the chain with zero magnetization has a different excitation structure. Let us emphasize that the fc-A string below 
half filling has a gap! Hence, the fc-A string does not influence the conformal dimensions [|j and does not lead to new 
low energy physics. The same statement holds true for all longer fc-A strings. 

(III) Next, we want to comment on eq. (7) of [Q. The authors of [0 do not require (7). They claim that this is the 
reason for "new" physics. All new physics they present is related to states containing one fc-A string (and several real 
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kj and real A-,,). The authors of [Q| do not argue the validity of their eqs. (5), (6) and (8). We show in appendix B that 
for one fc-A string (and several real kj and real A^) eq. (7) follows from (5), (6) and (8). This shows an inconsistency 
in the arguments of |^ and removes the basis for the modification of the physics of the one-dimensional Hubbard 
model. 

(IV) Finally, let us comment on the continuity argument on the bottom of the left column of the page 3 of [Q. 
The defect of the argument is the following: in the limit u oo [when the bound state (9) satisfy periodic boundary 
condition] the energy of the bound state is infinite, and it drops out of the Hilbert space. Continuity to finite u cannot 
be employed to count such states. 

We would like to thank H. Frahm for helpful communications. 



APPENDIX A: DERIVATION OF EQUATION (2) 

In this appendix we show that eqs. (20)-(22) in ^ are identical with the string Bethe equations of Takahashi and 
that (H) follows from these equations. 

Let us first cite the main results of : The rt-complex is parameterized by an n-string of the form 

'^S='/'i"^ + (" + l-2j)j , J = l,...,n. (Al) 

[Note the slight change of notation compared to [0 : We made the replacement (/)q"^ — *■ (and ^^""^ (fi^^^j ) which 
allows us to include scattering of two different n complexes of the same length, a enumerates the n-strings.] The 
S-matrix of an unbound particle with an n-complex is 

sinfc-0i"^ --i^ 



-,u(n) ^ 

sinfc — (fia^' +n\ 



= . , — (A2) 



and the S-matrix of an m-complex with an n-complex is 

_ ^1"^ - 4"^ ~\n~m\f ^\r^ - -jn + mYf U(^) „ _ ([^ - + 20f \ ' 

^i-'^i"' _ + In - m|f - + in + m)f t\ ^1"^ - -^^"^ + i\n - m\ + 21)^ ) ' 

Diagonalizing the transfer matrix leads to the following set of equations 

— 1 4 {n,a) 



~ As ~ ^ -r-r A-y — sm k 



//,A,-A, + f iiA,-sinfc, + f' 

AT" 

ei,-(.-).= Yl slfrJ^,l[Sl^^„ (A6) 

{m,b)^{n,a) j = l 

where 

q^"\<j)) = -2Rearcsin((/i + niu/i) + [n + 1)tt = -(arcsin((/) + niu/4) + arcsin(0 - niu/i)) + {n + 1)tt (A7) 

and 



Eqs (|A4| )-(|A6| ) are the "new" equations of Braak and Andrei ((20)-(22) or (BAG) in [|l|]). In order to show that 



(A4)-(A6) agree with Takahashi's string Bethe equations let us now adjust notations. Let 

, AT = 0i") , a = a , (3 = b. (A9) 
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Let us further introduce the functions g] 

X + i 



e{x) 



X — 1 

e 



\n — m\ 



X 



e - e- - 



|n — m| + 2 

X 



n + m — 2 



X 



n + m 



for n m, 



2n~2 



— I for n = m. 
2n) 



(AlO) 
(All) 



Then eqs. (A4)-(A6) turn into 



A-y — sin /cj 



(5=1 
A/" 



(5=1 



sin fcj — A5 

c7 



— A5 
2U 



(n.a) 



sinfcj - A'^ 
nf/ 



exp{— iL[arcsin(A'^ + n\U) + arcsin(A'^ — mU) + (n + l)7r]} 

A „ - sm fci \ T-r 

11 ^" 



sin few 



nU 



A 



(A12) 
(A13) 

(A14) 



The spin rapidities A-y in (A12) and (A13) may generally be complex. In order to obtain a set of equations which 
contains only real unknowns and which transforms into a set of linear integral equations in the thermodynamic limit 
we have to employ Takahashi's string hypothesis for A strings: As the number N of electrons becomes large the spin 
rapidities are driven to string positions characterized by their length n and their real center AJJ . Following Takahashi 
we will use the notation AJJ'-' instead of A^. AJ^'-' is the j-th spin rapidity involved in an n-A string with center AJJ, 



A'^'^' = A'^ + (n + 1 - 2j)iC/. 



(A15) 



Following again Takahashi let us assume that in the thermodynamic limit all A-, are grouped into strings with an 
accuracy of C'(exp(— (5A^)), where b is some positive number. Then eqs. (A12)-(A14) lead to 



\k^ L 



(n,a) 



AT" 



A" — sin kj 
nU 



(771,13) 



sinfc,--AS 
nU 

A" - A 



(n.a) 



sinfcj - A';^ 
nU 







u 



exp{— iL[arcsin(A'^ + iiiU) + arcsin(A'^ — niU) + {n + 



i=i 



nU 



(771,13) 



A'" _ A'™ 

a f3 

u 



(A16) 
(A17) 

(A18) 



Taking logarithms we arrive at the following form of the string Bethe equations, which is suitable for considering 
the thermodynamic limit, 



N-2M' 

E ' 



AJJ - sin k 

nU 



00 M„ 

kjL = 2^1, 

00 Mm 

= 2ttj^ + ^ ^ e„„ 

m=l /3=1 
N-2M' 



sin kj — A"^ 
nU 



A/', 



EE' 

n—l a—1 



sin kj — A' 



nU 



A" — A™ 

u 



L[arcsin(A'^ + inU) + arcsin(A';^ - mU)] = 2t:J''2 + E ^ 



A'^ — sin kj 
nU 



00 A(f;, 

E E 

m=l 13=1 



u 



(A19) 
(A20) 

(A21) 
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Here we assumed L = 2 x odd to be even. Ij, J", and J'" are integer or half-odd integer numbers, N is the total 
number of electrons, M' ~ X^^i'^-^rn ^(^) — 2arctan(a;), and 



\n — m\ 



29 



26 



20 



(9 



\n-m\ + 2 
■■ + 29 



20 



X 



2n - 2 



n + m — 2 
if n : 



X 



n + m 



if n ^ m, 



(A22) 



The branch of arcsin(.T) in (A21) is fixed as — 7r/2 < Re(arcsin(x)) < 7r/2. Af„ and M'j are the numbers of A strings 
of length n, and A' strings of length n in a specific solution of the system ( A19)-( |A21 ). The integer (half-odd integer) 
numbers in (A19)-(A21) have ranges 



L- 1 



L- 1 



N - 2M' 



1 / °^ 
Z\<^[l-N + 2M'^Y.^ 

^ \ m=l 



(A23) 
(A24) 

(A25) 



where tmn — 2 min(TO, rt) — (5mri- Each set of numbers {/j}, {>^'q} is in one-to-one correspondence with a set of 

rapidities {fc^}, {A^}, {A'^}, which in turn unambiguously specifies one Bethe eigenstate of the Hubbard Hamiltonian. 
Thus the ground state and all excited states can be constructed by specifying a set of numbers {/j}, {^q}, and 
then taking the thermodynamic limit. 

It is our aim to prove eq. (||) for the charge singlet excitation o ver th e half -fille d ground state. The ground state is 
characterized by TV = L, M\ = L/2. In this case the inequalities (A23) and (A24) lead to unique distributions of the 
quantum numbers Ij and J^, 



-(i + l)/2+j , J 



-(L + 2)/4-fa 



.,L/2. 



Eqs. (|A19| ) and (|A20|) reduce to 



Lfc, — 27r/i 



A„ 



L/2 

Q = l 

sin fc, 



U 



sin kj — Aq 
U 

L/2 

= 2^^a+E' 

13=1 



Ag - A^ 

2U 



a 



l,...,L/2. 



In the thermodynamic limit eqs. (A27) and (A28) turn into the well known integral equations |^ 

U 



11 f 

p{k) = - — I — cosfc / 

2Tr TT J _ 

a{K) = 



dA 



1 

2^ 



dk 



C/2 + (sinfc- A)2 
U 1 

TT 



C/2 + (sin k- A) 



a{A) 
dA' 



2U 



4t/2 + (A - A') 



(A26) 

(A27) 
(A28) 

(A29) 
(A30) 



for the densities p{kj) = \/{L{kj+i — kj)) and cr(Aa) = l/(L(Aa+i — Aq,)). 

The charge singlet is characterized by TV = L, Mi = L/2 — 1 and M{ = 1 ||,^. Thus M' = 1, and the number of 
unbound electrons is L — 2M' = L — 2. We will denote quantities which describe the charge singlet by a tilde. Eqs. 
( A24| ) and ( A25| ) uniquely determine the set {J^} and the number J' , 



J' 



-L/A + a a = 1, . . . ,L/2 - 1 



0. 



(A31) 



The set {Ij}, however, is not uniquely determined by the inequality (A25). There are (2) inequivalent such sets, 
which are parameterized by two vacancies I'l and I2 in the distribution (A26) of the numbers Ij, wh ich characterize 
the ground state. These two vacancies determine two charge rapidities ki and ^2 via eqs. ( |A27 ) and (A28). In 
the thermodynamic limit the charge rapidities densely fill the interval [— tt, tt], and ki and fcj become the two free 
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parameters of the charge singlet excitation. We see aheady at this stage, that there can not be a third free parameter, 
since J' = is fixed. For the charge singlet excitation eqs. ( A19 )-( A21 ) reduce to 

L/2-1 



Lk,=2^i,- J2 



sin kj — Aq 



a = l 



U 



sin kj — A' 
U 



, .7-l,...,L-2, 



L-2 

E' 



Aq — sin kj 



L/2-1 



Aa - Af 



/3=1 



2U 



, a = 1,..., L/2-1, 



L-2 



L(arcsin(A' + iU) + arcsin(A' - iC/)) = ^ i 



A' — sin kj 



(A32) 
(A33) 
(A34) 



Let us subtract (|A27D from ( |A32| ) for Ii ^ Ij, j ^ I, . . . , L - 2, and ( |A28| ) from ( |A33| ) for a = 1, . . . , L/2 - 1. Then, 
taking the thermodynamic limit, we obtain a pair of integral equations for the shift functions 



A 



Q + l 



A„ 



These integral equations are 



•^■^^ ^ 2 27r V ?7 



dA 



(=1 ^ ^ ■'^ 



dA' 



C/2 + (sinfc-A)2 
2C/ 



4C/2 



(A35) 

(A36) 
(A37) 



Here we have supplied an index "CS" to emphasize that we are dealing with the charge singlet excitation. The solution 
F^g{k) of (A36) can be found on page 526 of It is a single valued function of sinfc. In the derivation of (A3C) 
and (A37) we made use of the following elementary lemma. 



dk /(sin k) cos fc = 0, 

-7T 



(A38) 



which holds for arbitrary single valued functions f{x). This lemma follows from the identities sin(7r — fc) — sinfc and 

cos(7r — fc) = — cos fc. 

Note that there are still three free p arameters, fcf , k!^ and A', in (A36) and (A37). A' becomes fixed as a function 
of ki and k!^ by considering eq. (A34), which in the thermodynamic limit turns into 



L(arcsin(A' + ill) + arcsin(A' - iU)) 
'A' - sinfc 



L 



dk 9 



U 



— I dke 

27r 



A' 



U 



p{k) - 2 



dk 



U 



C/2 + (A' - sin fc)' 



;F^sik) cos k-J2(^ 



A'-sinfcf 
U 



A' - sin fcf 
U 



(A39) 



In oder to obtain the second of the above equalities we have used ( A29| ), (A36) and the lemma ( A3^ ). Let us calculate 
the integral 



ZTT 



dk 9 



A' - sinfc 
U 



on the right hand side of (A39). First note that /(O) = 0. The derivative of /(A') can be represented as 



/'(A') 



1 



dk 



U 



C/2 + (A' - sinfc)2 
where the contour of integration is the unit circle. Let 



Re 



27ri 



dz 



z2 + 2ziU - iA') - 1 



(A40) 



(A41) 
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p{z) = + 2z{U - iA') - 1 = {z - zi){z - Z2). 



(A42) 



We see that the poles of the integrand in ( A41 ) are related as 

Zi = -1/Z2. 



(A43) 



Thus, only one of these poles, say Zi, lies inside the unit circle. Using (A43) we obtain I' {A.') as a function of zi, 

2 



/'(A') 2Re 



Zl + z^ 



(A44) 



Let us parameterize the pole at zi as zi — e'". Since zi is located inside the unit circle, Im(a) > 0. Using p{zi) — 
we find that 



sin a = A' + iU , lm{a) > 0. 
This equation fixes a modulo 2tt. Now U > and Im(a) > 0, and therefore 

U 



cosRe(a) 



sinh Im(a) 



> 0. 



(A45) 



(A46) 



We conclude that — 7r/2 < Re(Q:) < n/2. Thus (see the definition below (A22)) a = arcsin(A' + i[/). Integrating (A44) 
with respect to A' and using 7(0) = to fix the integration constant we arrive at 



/(A') = 2Re(Q;) = 2Re(arcsin(A' + it/)) = arcsin(A' + iU) + arcsin(A' - it/). 



We may now insert this result into eq. ( |A39D . This yields 



A' -sinfcj' 



A' - sin k'^ 



U J \ u 
Dividing by 2 and taking tan gives 

2 A' = sinfc^' + sinfc2. 
So we have accomplished our task to show that (2) follows from eqs. (20)- (22) of 



(A47) 



(A48) 



(A49) 



APPENDIX B: 



Let us show that for one k-K string (and several real kj and real A^) eq. (7) of j^] follows from eqs. (5), (6), (8) 
of In case of a single fc-A-two string the numbers iV" and Af" of are — N ~ 2, A/" = A/ — 1, where N 
is the total number of electrons and M is the total number of down spins. There are TV — 2 real charge rapidities 
ki, . . . , fcAr_2, which correspond to unbound particles. The string is characterized by two charge rapidities A:"*" and k~ 
and by one spin rapidity A. k^ and A satisfy 



sin A;+ = A - iit/4 + e+ 
sin/c^ = A + iu/4 — e'~ 

where e+ and become small as the length L of the lattice becomes large. Eqs. (5)- (8) of 0| are 



\~f — Xs — iu/2 
Xj — Xs + iu/2 

M-l AN • /r, 

nA — A5 — m/2 
K-Xs+ iu/2 



A — sin fc, 



iu/4 -j-j- Xs—smkj 



iu/4 



M-l 

n 



A — sin kj + iu/4 ^ Xs — sin kj + iu/4 

|^ A,-sinfe, i./4 

A-, — smKo + m/4 
.7=1 ' 



,iV-2, 



-n- A — smfcj — m/4 



J{k++k-)L 



-n 

A/-1 

n 



e A — sin fc," + iM/4' 

7 — 1 ' 

Aa - A - iu/2 



(5=1 



Xs- A + iu/2' 



(Bl) 
(B2) 



(B3) 
(B4) 

(B5) 
(B6) 
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The ratio 



£+ A — sin fc+ — iu/4 
e~ A — sin + iii/4 



(B7) 



is denoted by e"^ in [||. We want to show that ( |B5| ) algebraically follows from (B3), (B4) and (B6). First note that 
momentum conservation implies that 



exp (iL (Ejlf A:, + fc+ + fc-)) = 1. 



Multiplying all eqs. ( |B4|) we get 

7V-2 A/-1 



. ^ \\ Ai-sin%+iu/4 " 1^ ^ Ai - A^ + m/2 " ' 

7^5 



Now let us multiply all eqs. (B3). Taking into account (B£) we get 

= TT A - sin % -1^/4 
A - sin kj + iu/4 ' 

7 — 1 J I 



Let us multiply this by (B6) and use (BS). Then 
1 = exp (iL (EjI^i' fcj + /c+ + fc 



-p-i- A5 — A — m/2 -p-i- A — sm kj — lu/i 
~ Xs -A + m/2 A - sinfc,- +iM/4' 

(5=1 ' 1 = 1 ' 



(B8) 



(B9) 



(BIO) 



(Bll) 



which is equivalent to (B5) (or (7) in [hj). 
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